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Abstract. Driving nanomagnets by spin-polarized currents offers exciting prospects in 
magnetoelectronics, but the response of the magnets to such currents remains poorly un- 
derstood. We show that an averaged equation describing the diffusion of energy on a graph 
captures the low-damping dynamics of these systems. From this equation we obtain the 
bifurcation diagram of the magnets, including the critical currents to induce stable preces- 
sional states and magnetization switching, as well as the mean times of thermally assisted 
magnetization reversal in situations where the standard reaction rate theory of Kramers is 
no longer valid. These results agree with experimental observations and give a theoretical 
basis for a Neel-Brown-type formula with an effective energy barrier for the reversal times. 



Manipulating thin-film magnetic elements with spin-polarized currents besides external 
magnetic fields [I2| has generated a lot of recent interest in applications to magnetoelec- 
tronic devices that offer low power memory storage without the use of moving parts P]. 
Understanding the response of the magnet to such currents is nontrivial, however, because 
they apply a nonconservative force, called spin-transfer torque (STT), on the system. Like 
other nongradient systems with no Lyapunov function, the phase portrait of nanomagnets 
in the presence of STT can be quite complex, and include limit cycles or chaotic trajecto- 
ries besides fixed points. When the applied fields and/or currents are nonstationary, or in 
the presence of thermal noise, the situation is even worse. In particular, Kramers' reaction 
rate theory [131 ITO] is no longer applicable and the Neel-Brown formula pj for the mean 
magnetization reversal time is not valid since there is no well-defined energy associated with 
STT. 

Nanomagnets typically operate in a regime where the nonconservative parts of the dynam- 
ics, including the effects of damping, STT, and thermal noise, act on time-scales that are 
much longer than that of the energy-conserving Hamiltonian part. Trajectories remain close 
to periodic Hamiltonian orbits for a long time, and slowly drift from one orbit to another 
due to damping, STT, and thermal noise. This separation of time scales can be exploited, 
using averaging techniques developed by Freidlin and Wentzell [H [8] (see also [201 S!); to 
reduce the dynamics to that of an energy diffusing on a graph. We show here that this 
reduced description permits to explain the features of nanomagnets subject to STT that are 
observed experimentally. Specifically, we obtain the full bifurcation diagram of the system at 
zero temperature and determine the critical spin-polarized currents needed to induce stable 
precessional states [15], |3| and magnetization switching [151 [19]. At finite temperature, we 
calculate the mean times of thermally assisted magnetization reversals [151 ^S\, and give 
expressions for the effective energy barriers conjectured to exist [151 [HI [D HSl [II]. These 
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results are complementary to those obtained in [6|, using the geometrical Minimum Action 
Method (gMAM) [11], for situations with small thermal noise and stronger damping. 

We will focus on magnetic systems in which the magnetization has constant strength Ms 
in the direction of a unit vector m{t) = {mx(t),my(t),mz{t)) whose evolution is governed 



This is the standard stochastic Landau-Lifshitz-Gilbert (LLG) equation written in non- 
dimensional form (details of the nondimensionalization can be found in the Appendix A), 
with an additional STT term, m x (m x ajrUp), modeling the transfer of angular momentum 
to the magnetization from the electron spin in a polarized current of strength aj directed 
along unit vector nip [21j. For simplicity, here we take rUp = (1,0,0) and a constant 
strength aj, but these could straightforwardly be generalized to any direction and a time- 
varying strength. The other terms in ([1]) are standard: a is the damping parameter and 



hcff = —VmE + y Y^'ni't) is the effective field, which in turn is the sum of the negative 

gradient of the energy per volume, E{m), and a term accounting for thermal effects with r}{t) 
being three-dimensional white-noise and e = ksT / iiqM'^u the non-dimensional temperature 
(z/ is the magnet volume and /zq the permeability of free-space). Here we take E{m) = 
f3ymy + Pzml — hxTUx with (3y < (3z = 1/2, corresponding to biaxial anisotropy with a planar 
applied field of amplitude hx in the x-direction. 

The presence of the STT term in ([T]), which is nongradient and nonconservative, compli- 
cates the analysis of this equation even in the absence of thermal noise (e = 0). In particular, 
the magnetic energy E{m) is not a Lyapunov function for the system, and it is unclear if 
the STT term provides additional damping, driving, or something entirely different. To 
understand the effect of this term, we take advantage of the separation of time scales that 
arises when both the damping and the strength of the polarized current are weak, a -C 1 
and aj ^ 1. In this regime, m moves rapidly along the energy conserving Hamiltonian 
orbits in Fig. [I](a) and drifts slowly in the direction perpendicular to these orbits. This slow 
motion can be captured by tracking the evolution of the energy E{m) along with an index to 
distinguish between disconnected orbits with the same energy. This information is encoded 
in the graph shown in Fig. [I](b), whose topology is directly related to the energy function, 
E = /3ymy + Pzml — hxirix, and changes based on its form and values of parameters. For 
example, when l/i^l < 2/3^ the graph has four branches, as shown in Fig. [D^b), which meet 
at the saddle point of the energy that corresponds to the homoclinic orbits connecting the 
two green points on the surface of the sphere in Fig. [l](a). We will use the indexes 1 and 2 
(3 and 4) for the lower (higher) energy branches in Fig. [I](b), which correspond to orbits on 
the front-right and back-left (top and bottom) of the sphere in Fig. [l](a), respectively. 

To deduce the effective dynamics on the graph when a and aj are small, we follow Freidlin 
and Wentzell [9] to remove the direct dependence on m{t) from E = VmE ■ rh by time- 
averaging the coefficients in this equation over the Hamiltonian orbits of constant energy 
depicted in Fig. [l](a). The resulting averaged equation for the energy on branch j is (see 
Appx. B) 



by 



(1) 



m = —m X hcff -|- m X (m x (— ahefr -|- ajnip)) . 




(2) 



E 
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Figure 1. (color online) (a): Hamiltonian orbits (blue lines) of the unit mag- 
netization vector solution of m = —m x VmE along with fixed points of this 
equation that are also critical points of the energy: minima (red dots), saddle 
points with energy Eb (green dots), and maxima (blue dots), (b) and (c): 
Graphs (not drawn to scale) in which each energy point along the edges corre- 
sponds to an orbit of m shown in (a) with this energy. The numbers indicate 
the label for each energy branch. Branch 1: rrix > —hx/2Py and E < Et, 
branch 2: rrix < —hx/2Py and E < Eb, branch 3: rUz > and E > Ei, and 
branch 4: < and E > Eb- The ends of each branch correspond to the 
fixed points in (a) and the circles indicate the location and stability of the fixed 
points of E = —aAj{E) + ajBj{E): filled circles are stable and open circles 
are unstable. The graph in (b) corresponds to a situation where the energy 
minima are stable fixed points. The graph in (c) is a representative case when 
a solution of exists, leading to a new stable fixed point ai E = Eq] this 
new fixed point corresponds to a stable limit cycle like the one shown in red 
in (a). 



written in Ito's form, where ^(t) is a ID white-noise and 

ME) = A{Pl (ml)^ + Pi {ml)^ - E') - AEh. (m.)^. + hl{l - {ml) ^ ) 
(3) B,{E) = 2E{mx)^ + K{l + {ml)^) 

C,iE) = Py + P,-3E-2K{mx)^ . 

Here denotes the time-average over one period along the orbit with constant energy E 
corresponding to branch j of the graph in Fig. [T](b). As we show in the Appx. C, the averages 
in (|3]) can be evaluated asymptotically near the critical points, and this information turns out 
to be sufficient to calculate the bifurcation diagram and the mean times of magnetization 
reversal that we obtain below. Away from the critical points, the averages must be 
evaluated numerically, which we do by using a symplectic implicit mid-point integrator to 
evolve m via rh = m x VmE along an orbit with prescribed energy to compute the time 
averages. Note also that ([2]) requires a matching condition where the branches on the energy 
graph meet |9]; these conditions are discussed in the Appx. D. 

Next we use the reduced equation ([2]) to obtain the bifurcation diagram of the system at 
zero temperature, e = 0, and determine the fixed points of E = —aAj{E) + ajBj{E) and 
their stability. The coefficients Aj{E) and Bj{E) encode the separate effects of the damping 
and the STT on the energy, respectively, and it can be checked that they are both zero at 
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the critical points of the Hamiltonian. The energies at these points are 

Ea,i = -K and Ea,2 = 

corresponding to the two energy minima on the lower branches 1 and 2 where m = (±1, 0, 0), 
respectively; 

corresponding to the saddle point in energy where all four branches meet where m = 



(-V2/3„±^l-/i2/4/32,0); and 

corresponding to the two energy maxima on the upper branches 3 and 4 where m = 
{-K/2I3,, 0, ±^1 - hi/ API). These critical points can merge and disappear when the ap- 
plied field crosses the critical values hx = ±2/3^ and hx = ±2/32. In addition, only the two 
energy minima can ever be stable, and one of them looses stability when another nontrivial 
fixed point in energy, Eq, appears on one of the energy branches, where the energy lost by 
damping, —aAj^Eo), is exactly compensated by the energy gained by STT, ajBj^Eo): 

(4) -«A,(Eo) + aj5,(Eo) = 0. 

The stable fixed point at Eq does not corresponds to a stable fixed point of the original 
dynamics, but rather to a stable limit cycle (precessional state), see Fig. [T]for a schematic 
illustration. The location and stability of the fixed points identified above are shown in 
Fig. |2]^a) as a function of aj for a fixed value of hx, highlighting that magnetization reversal 
can be achieved by varying the strength of the spin-polarized current, aj. 

We can also calculate the full bifurcation diagram shown in Fig. |2]^b), which is remarkably 
similar to the experimental one (see Fig. 2a in [H]). One of the energy minima looses its 
stability and the precessional state appears when Ea i or Ea 2 solves (jl]), i.e. when aj is given 
by(j = l,2) 

A(x) 

(5) aj = a lim ^ = aaj{(5y + (3^ + (Tjhx), 

where cJi = 1 and a2 = —1. The corresponding boundaries on the bifurcation diagram are 
shown as dashed lines in Fig.HJ^b). The limit in was obtained using asymptotic expansions 
of the coefficients; details can be found in the Appx. C.l. The precessional state exists in 
the region between the dashed and the solid lines in Fig. [2]^b). Beyond these solid lines only 
one stable state remains. This occurs when Ef, solves dl]), meaning that the strength of the 
current required to induce switching is (j = 1,2) 



aj = a lim , , = A 

^3 



' a,h,hx + T^^,Py{l-hl/APl) 



where h, = {^S/K + ' Py)hx/2f3yf3x)'', d, = 1 - {^S/K - cTjhx/2f3,)\ and S = 

{Pz — Py){Pz — hl/4:l3y); details can be found in the Appx. C.2. 

Next we study thermally induced magnetization reversal. To this end we use the reduced 
system in ([2]) with e > to calculate the mean transition times (i.e. dwell times) between the 
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Figure 2. (color online) (a): Bifurcation diagram as a function of aj at the 
fixed value of = 80.64 G: the stable (solid line) and unstable (dashed line) 
fixed points of E = —aAj{E) + ajBj{E) are shown for both energy branch 
1 (blue) and 2 (red). (The remaining two unstable fixed point with higher 
energy are not shown.) (b): Bifurcation diagram as a function of aj and 
hx (converted to units of Gauss by fioMslO'^, Mg = 954930 A/m, to easily 
compare with experimental data shown in Fig. 2a in [15j). The dashed lines 
on either side of region S correspond to the current required to first initiate a 
stable precessional state (aj given by ([5])) while the solid lines correspond to 
the current required to induce switching (aj given by ([6])). Beyond the solid 
lines, only one stable fixed point remains: in region Tp (T^) it is m^, = 
(m^; = — 1) and for beyond the region shown, only one lower energy branch 
remains. The stars indicate where the mean thermally induced switching times 
from branch 1 to 2 and branch 2 to 1, computed via ([7]), are equal. For both 
figures, the values a = 0.003, (3y = 0.0654, and (3z = 0.5 are used. 



stable fixed points of the deterministic dynamics identified before. The mean time Tj{x) to 
transition from energy x on branch j = 1, 2 to the fixed point on the other branch satisfies 

[ — aAj{x) + ajBj{x)+2aeCj{x)~\Tj{x) 

+aeAj{x)T^'{x) = -1 

with a matching condition (see Appx. D) to prescribe transitions through the center node 
of the graphs shown in Figs. [U^b), (c) as well as an absorbing condition at the target state. 
Equation ([7]) is valid at any temperature and its solution can be expressed in terms of 
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integrals involving the coefficients Aj{x), etc. Evaluating these integrals numerically leads 
to the results shown in Fig. [51 We can also evaluate these integrals asymptotically in the 
limit when the temperature is small (e <^ 1), in which case they are dominated by the 
known behavior of the coefficients near the critical points. These calculations are tedious 
but straightforward and reported in the Appx. E. In situations where the system transits 
from the stable minimum Ea^i or Ea^2 on branch 1 or 2 to the stable point (minimum Ea^2 
or Ea^i or precessional state Eq) on the other branch we obtain {j = 1, 2) 

(i-^)(i?,-K,,)A 

(8) r„- ~ 



whereas in situations when switching occurs from the stable precessional state Eq we obtain 
(j = 1, 2 depending on whether Eq is on branch 1 or 2) 

il-^)iE,-Eo)/e 

(9) r - 

Here A^, defined in ([S]), is the critical current to induce switching at zero temperature, and 

7i = 4d2(/3.'-/3j)-&2(4/3j + /^^) 

with bj and dj defined after The results in and agree with the experimental 
observations [151 IS] and the theoretical predictions [H that the effect of STT on the 
dwell times can be captured via a Neel-Brown-type formula with an effective energy scaling 
linearly with the current strength. We stress, however, that these previous theoretical works 
had to assume the existence of such a formula, whereas ([8]) and (Q fall out naturally from 
the asymptotic analysis, and give explicit expressions not only for the effective energy but 
also the prefactors and their dependency on the strength of the current producing STT. 

In summary, we have shown how the dynamical behavior of nanomagnets driven by spin- 
polarized currents can be understood in the low-damping regime by mapping their evolution 
to the diffusion of an energy on a graph. We thereby obtained the full bifurcation diagram of 
the magnet at zero-temperature as well as the mean times of thermally assisted magnetization 
reversal. These results agree with experimental observations and give explicit expressions 
for the dwell times in terms of a Neel-Brown-type formula with an effective energy, thereby 
settling the issue of the existence of such a formula. We carried the analysis for micromagnets 
that are of the specific type considered by Li and Zhang [16], but the method presented in 
this paper is general and can be applied to other situations with different geometry, applied 
fields that are time-dependent or not, etc. It can also be applied to systems in which the 
magnetization varies spatially in the sample. In these situations, the graph of the energy 
will be more complicated, but the general procedure to reduce the dynamics to a diffusion 
on this graph remains the same. Such a study will be the object of a future publication. 

The research of E. V.-E. was supported in part by NSF grant DMS07-08140 and ONR 
grant N00014-11-1-0345. 
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Figure 3. (color online) Contour plot of the mean first passage time (con- 
verted to seconds by 7/ioM^/(l + a^^, 7 = 2.21 x 10^ m/A-s, M, = 954930 
A/m) from the fixed point on energy branch j = 1, to the fixed point on branch 
2. (The times from branch 2 to 1 would be the same figure, rotated 180°.) 
These times are computed using ([7]) and plotted as a function of (converted 
to units of Gauss by ^^MglQ'^) and aj. The shape of the plotted region is 
identical to that shown in Fig. |2]^b). Between the dashed line (aj in ([5])) and 
solid black line (aj in ([6])), the starting point is Eq (Eq. (jl])); everywhere else 
it is the minimum energy, -Ea,i- The values a = 0.003, e = 0.01, (3y = 0.654, 
and = 0.5 are used. Switching times greater than 25 years are all colored 
in dark red. 



Appendix A. Nondimensionalization of the Governing Equation 

We focus on magnetic systems of the type considered in Li and Zhang ^16j. They consist 
of a thicker fixed-magnetization layer, separated from a thinner layer whose magnetization 
we wish to model. This free-layer's magnetization vector M = (M^., M^, M^), assumed to 
have constant strength |M| = Mg {Mg is the saturation magnetization), is governed by the 
(dimensional) stochastic Landau-Lifshitz-Gilbert (LLG) equation 

dM 'y*a 
(10) = -7*M X Heff - -^M X (M X H,ff) 

at Ms 

where 7* = 7/(1 -|- a^), 7 is the gyromagnetic ratio (units rad/s-T) and a is the non- 
dimensional Gilbert damping constant. The effective magnetic field (units T), 



(11) Heff=-VM^+ J^^r7(t), 

is the sum of the 3-component white noise, T]{t), and the gradient of E, the energy per 
unit volume (units N-m/m^). The noise amplitude, derived from the fluctuation-dissipation 
theorem, will be discussed later, u is the volume of the magnetic element, its temperature 
is T, and the Boltzmann constant is ks = 1.38 x 10~^^N-m/K. The energy per unit volume 
varies across different materials. Here, we use 

(12) E = /.oH.., ■ M + ^M^^ + ^Ml 
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the sum of the contributions from the external field, Hext, and the anisotropy in the magnet 
along the y direction, given by Hj^ (hard axis z, easy axis x). The permeability of free space 
is /io = 47r X 10"^ T-m/A. 

In order to non-dimensionalize Eq. (ITU]) , we first write M = M^m where m = {m^j., niy, rriz) 
is a unit vector, E = ^qM^E where E is the non-dimensional energy per volume, and 
Hcff = fioMgheff. With these substitutions we obtain 

= -YfioMsin X hcff - Yf^oMsam x (m x h^s) 

at 



(13) hcff = ~VmE + xl^^^^rjii) 

'H-ext 1 -f^A: 2,1 2 

After non-dimensionalizing time to t = t/'~f*fiQMs, the resulting non-dimensional equation is 

U /UN 

-m X hcff — am x (m x heff) 

UL 

^^^'^ ' 2akBT 

r,it). 




(1 + a^)i2oM^i^ 



Note that the appearance of (1 + a^) was because 7, not 7* appeared in the denominator of 
the white-noise amplitude factor. 

To Eq. (fTi|) . we add the spin-transfer torque (STT) term, ajm x (m x nip), describing 
the response of the system to a spin-polarized current with (non-dimensional) strength aj 
generated from the fixed magnetic layer, with magnetization unit vector nip. The non- 
dimensional strength aj = 7*/io?7(6')/iB//ei/, contains the dependence on the current, J, and 
the structure of the fixed and free magnetic layers through 7]{6) = q/{A + B cos{6)) where 
cos(6') = m • nip [2T]. The current strength is divided by the volume, u, of the free layer, 
so as to produce a force per volume, matching the energy per volume contribution already 
in Eq. f|T^ . We take the value of aj to be constant here, though our results could easily be 
extended to the time-varying case. 

For simplicity, we consider both the spin-polarized current direction, nip, and the external 
field, Hext, to be directed along the x-axis: nip = (1, 0, 0) and iiext = {hxMg, 0, 0). The final 
non-dimensional LLG equation that appears as Eq. (1) in the text is 

(15) ^ = -m X hcff - m X (m X (ahcff - ajnip)) 



where heff = -V^E + sj ^vit) with e = ksT/f^oM^u and 
(16) E = j3yml + (5^ml - h^m^ 

and where fiy = Hk/2Ms < Pz = 1/2. Recall that Ms is the constant amplitude of the 
dimensional magnetization vector and u the volume of the magnetic element. 

We return to consider the amplitude of the noise, ^2ae/{l + a^) in Eq. taken 
from the paper by Kohn, Reznikoff and Vanden-Eijnden [13]. The parameter e is the non- 
dimensional temperature; it is the thermal energy, ksT, non-dimensionalized by noM^u, 
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where v is the volume of the magnet and fioM^ is the energy per volume scaling. The a 
appears in the noise term so that the equilibrium distribution is the Gibbs distribution, 

p(m) = Z-le~^('")/^ 

where is the normalization factor; it is the integral of e~^^"^^^^ over the surface of the 
m-sphere. This fact is shown in Appendix A of the paper by Kohn, Reznikoff and Vanden- 
Eijnden [13]. It was done by switching to the {6, z) representation of the magnetic vector to 
represent this two degrees of freedom system by two variables. From the non-dimensional 
form of the noise term, we worked backwards through the non-dimensionalization procedure 
to obtain the dimensional noise amplitude in Eq. ( ITOj) . 



A.l. Converting to Ito form. We convert the Strotonovich SDK flTB]) into an Ito SDE, 
which will be needed in the next section to determine the averaged evolution of the energy. 
First, we write the Strotonovich SDE in the form 



„x dm / 2ae „ , , 

(17) -— = He - aarf + ajHp + W — — 7:Br]{t) 

dt V 1 + a 



where the conservative term is 



2(/32 - f3y)mymz 
-m^(2/3^m^ + 



the damping term is 



- ml) - 2m^{l3yml + I3zml) 
= I my{h^m^ + 2/3y(l - mf) - 2/3^m2) 
mz{h^m^ - 2(3yml + 2/3^(1 - mfj) 

the spin-torque transfer term is 

ml-1 



S.p — I TflxTfly 



and the diffusion matrix is 




a;(l — m^) — arrixmy —my — am^m^, 

B = I — — am^my a{l — my) 

my — am^m^ —m^ — amym^ 

In order to convert to the Ito form, the drift term obtains the correction 

a/ = — ^^-'S^ Bh-jduBi^ = —2aem, 
l + a^2^ ^ ^ 

making the Ito SDE for the magnetization direction 



dufx / 2ctc 

— - = a^ - aaa + a jap - 2aem + W — — :^B'n{t). 
dt V 1 + 
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Appendix B. Energy Equation 

Here, we derive the energy equation (2) in the text. Using the rules of Ito calculus, we 
compute E = VmE ■ rfi to obtain the Ito SDK for the evolution of the energy in the form 

(19) E = -aA{m) + ajB{m) + 2atC{m) + V2^t^ A{m)i{t) , 

where ^{t) is ID white noise, the strength of which is computed from the combination 
of the three strengths of the three independent components of the white noise, ^/(t), in 
Eq. f|T5l) . Notice since the term in Eq. f|T8|) conserves energy, it has no corresponding 
term in Eq. f|T9l) . The remaining terms in Eq. f|T8l) have corresponding terms in Eq. f|T9|) : 
the dissipative term, a^, leads to 

A{m) = A{Plml + _ _ ^Eh^m^ + hl{l - ml), 

the spin-torque transfer terms, a^, leads to 

B{m) = 2Em^ + h^{l + ml), 

and the correction term for Ito calculus {dj indicates partial derivative with respect to the 
j'th element of m), 

J2[SB%d,d,E{m) = 2ae(/3,(l - mj) + - m^), 

together with the contribution from a/ give 

C{m) = /3y + /3,-3E- 2Km,. 

The three separate contributions to the noise term, 

b^{m) = 2{(5z - (5y)mym^ - a[h^{l - ml) + 2m^{(3yml + /^^m^)] 
hy{m) = -m^{h^ + 2f3^m^) + a[hr,m^my + 2my{(3y{l - m^) - (3^ml)] 
b^{m) = Kmy + 2(3ym^my + a[hr,m^m^ + 2m^{f3^{l - ml) - fSymy)], 

simplify so that ^/^^^^-T^TT^ = a/(1 + a'^)A{m). To remove the explicit dependence on the 
magnetization vector m from Eq. (fT9|) . following Freidlin and Wentzell [9], we average the 
coefficients appearing in the backwards Kolmogorov equation for SDE ( |T9l) over one period, 
Tj{E), at constant energy. 

The subscript j = 1,2,3,4 indicates that the average corresponds to one connected orbit of 
m with constant energy E on branch j of the energy graph (Fig. 1(b) in the text). The 
resulting averaged coefficient backwards Kolmogorov equation corresponds to the averaged 
coefficient SDE (2) in the text. 

Appendix C. Asymptotic Expansions of the Averaged Energy Coefficients 

While explicit expressions for the averaged coefficients are unknown, we work out asymp- 
totic expansions near the two minima and the saddle point of the energy. These expressions 
show how the coefficients Aj{E) and Bj{E) go to zero near these points, and are used to 
determine the accessibility of these critical points in energy. The expressions are also used 
to determine the critical spin polarized current aj for the emergence of a stable precessional 



DYNAMICS OF STT MAGNETS USING ENERGY DIFFUSION ON A GRAPH 11 

state (Eq. (5) in the text), and for the critical current to induce switching of the magne- 
tization vector (Eq. (6) in the text). In addition, they are used in Sec. [D] in determining 
the switching probabihties between the lower energy branches and in Sec. [E]in asymptoticly 
approximating the mean switching time. 

C.l. Approximation Near the Minima. In order to determine the scaling of the averaged 
coefficients Ai{E), Bi{E) and Ci{E), near the energy minimum on branch 1, mg = (1, 0, 0), 
we create a series expansion about this point for the solution to the Hamiltonian system, 

/ 2{/3^ - /3y)mym^ 

(20) m = — m X V-E" = —■mz{2(3zmx + hx) 

\ my{2l3ym^ + h^) 

then compute the averages exactly as a function of the energy. We consider the expansion 
for m, 

m = (1, 0, 0) + (5mi + S^m2 + 0(6^), 

(where the subscripts denote terms in the expansion, and not energy branches) and determine 
the averaged coefficients up to 0(5^), requiring the first non-zero term above the zeroth term 
for each component of m. This expansion must also satisfy |mp = 1 up to order S"^ 



m^l = (1 + 5mi,^ + 6^m2,x) + ^ « + ^ z + 0(6' 



,y 

= 1 + 25mi,^ + 5^{ml.^ + 2m2,^ + m\y + J + 0(5^). 
To satisfy this constraint at 0(5), we must have that vtix^ = 0, and then at 0(5^), 

(21) m2,x = ~{mly + ml^). 
The 0{6) expansion of the system (!20!) is 

/ \ 

(22) mi = -2/3,-/1, mi. 

\0 2f3y + h^ / 

This is consistent with mi^x = obtained from the length constraint on m above. The 
remaining two components are 

rrii^y = Ci cos cot 

rriiz = Ci — — sm ut 

where to = a/ {2(3y + h^.) (2/3, + h^) and the 0{5) initial condition is assumed to be of the 
form mi(0) = (0, Ci, 0). No higher order terms are needed for the y and z components, but 
we must consider the next order term for the x component. 

For the component m2,x at 0(5^), the expansion of system (120|) gives 

rn2,x = 2(/3, - /3y)mi,j^mi,, . 

The solution to this is consistent with the definition of m2^x in (l2T!l need to satisfy the length 
constraint. By taking the derivative of m2^x given by (1211) . 

m2,x = -mi^yrrii^y - "^i,,mi,,, 
and then substituting the 0{6) DEQ for rfii y and mi we see that 

m2,x = '2{/3z - (3y)mi^ymi^z- 
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Combining these results we have the expansions 



2 

my ~ 5ci cos (jjt 



cos ujt H — sm ut 

2/3, + K 



where u = a/ (2/3^ + hx){2[5z + h^). These solutions correspond to a trajectory with initial 
condition (1 — 5^c^/2, 5ci, 0) and constant energy 

(23) E = -hx + 5^cl{Py + hj2). 

To determine the averaged coefficients, we average the functions of m over one period, 
T = 2tt /u (note that T does not depend on the energy in this expansion), 

U{m)) = ^ j^^ f{m{t))dt. 
Taking Ci = 1 (without loss of generality), we have 







2Py + K 




213, + hx 






2(3y + hx 




2(3, + hx 



2V2/3. + /^., 

To write the averages in term of the energy, we solve for 5^ as a function of E from Eq. 

52 E + hx 



2 2Py + hx' 



and obtain 



(24) 



Ai{E) ^ 2{[3y + [3, + K){E + hx) 
B^{E) ^ 2{E + hx) 



C.(i.)^/3, + /3. + /..-(l + ^ + ^)(E + M. 

These approximations are in good agreement with the numerically obtained averaged coeffi- 
cients, shown in Fig. |H The limit in Eq. (5) in the text is obtained from the ratio of Ai{E) 
to Bi{E). 

We repeat the above procedure for energy branch 2, where the minimum is m = (—1, 0, 0), 
using the series expansion 

m = (-1, 0, 0) + 5mi + 5'^m2 + 0{5^). 
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Figure 4. Comparison of the averaged energy equation coefficients on energy 
branch 1 computed by averaging the numerically integrated trajectories given 
by rh = — m x VE (blue solid line) to the asymptotic approximations in 
Eq. fl24|) near the minimum energy (green dashed line) and the asymptotic 
approximations in Eq. (l36|l near the saddle point in energy (cyan dash-dot 
line) for non-zero h^- The values I3y = 0.06, = 0.5 and = 0.03 were used. 
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From the constraint on the length of m we similarly find that mi x = and 

The 0{S) expansion of the Hamiltonian system has the solution 

rrii^y = cosut 

"^M = -\/ifcfesinu;t , 



where u = \/{2(5y — hx){2/3z — hx) and we have taken an 0{5) initial condition (0,1,0). 
Combining and averaging as above, using 

E- hx 



2 2(3y - K 



we obtain 



(25) 



A2{E)^2{Py + P,-hx){E-h 
B2{E) ^ -2{E - hx) 



for the averaged energy coefficients on energy branch 2. The limit in Eq. (5) in the text is 
obtained from the ratio of A2{E) to B2{E). 

C.2. Approximation Near the Saddle Point. We repeat a similar process to the one in 
the above section, and determine the scaling of the period of the orbit, which goes to infinity 
as the energy approaches its saddle point value, by again using a series expansion of the 
solution. To determine the approximation of the averaged coefficients, we must also consider 
their value along the entire orbit in order to compute the time-average integral. The integral 
is dominated by the values the components of m take along the homoclinic orbit connecting 
the two fixed points. We consider this integral last. 

Orbits with either energy slightly less than or slightly more than the saddle point value, 
Eb, spend the majority of their time near the two fixed points. To approximate the period, 
we expand the solution of the Hamiltonian system, Eq. f l20l) . about one fixed point and 
determine the time it takes the trajectory to leave an 0(1) size box around this fixed point. 
We must handle the two cases, hx = and hx 7^ 0, separately, but in the end we obtain the 
same scaling of the period with the energy in both cases. 

First, we consider the simple case when hx = 0, and consider the expansion of the solution 
about one of the fixed points, mg = (0, 1, 0), in the form 

m = (0, 1, 0) + 5mi + S^m2 + 0(5^). 

From the constraint on the length of the vector, 

1 = |m|2 = 1 + 25mi,y + S^{2m2,y + + m\ y + J + 0(5^), 

we see that at 0[S) it must be that vrix^y = 0, and at 0(5^), 

(26) m2,y = -^{ml^ + ml,) . 

The constant energy of this expansion can be found in terms of the initial conditions, and is 

E = Py + 2Py6^m2,yiO) + Pz6^mlM + 0(5^). 
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Due to the constraint, we only need to find the 0{6) solutions, then construct m2,y from 
Eq. ([26]). 

The expansion of the Hamiltonian system (!20|) leads to the 0{6) set of equations, 

rhi y = 



which is solved by 



mi J. = cosh ut 



where u"^ = 4:(3y{f3z — (3y). We have used the initial conditions (5, 1 — |5^,0) and find the 
energy scales with 6 as, 

E = /3y-l3y6^ + 0{6''). 

The approximate solutions ( 1271) do not lead to complete orbits as in the previous section, 
rather, we can estimate the period by computing the time for the trajectory to leave a box 
around the saddle point. Specifically, we approximate the period as quadruple the time it 
takes the m^, component to increase to some 0(1) value c away from its initial value. (We 
would obtain the same scaling if we chose to use the component instead.) Note that we 
are exploiting the symmetry of the system: there are two saddle points which the trajectory 
spends equal time near, and the trajectory spends an equal amount of time entering the box 
and traveling to the initial condition as it does leaving the box. Therefore, we wish to solve 

T 



5 cosh I u- 
V 4 

for the period T, of the orbit, which we expect to be large. We approximate cosh.uT/4 
(^^o^T/i ^ e-'^^/^)/2 as e'^^/V2 and obtain 

(28) r = llog-l + 0(l) = ^log-l^ + 0(l), 



where u = 2^J [5y{(5z — l3y) and the saddle point energy. Eh = l3y. 

We repeat the above procedure for the second case, hx 7^ 0. Again, we expand the solution 
about the fixed point, (—a, \^l--a?, 0) where a = hx/2l3y, in the form 



m = (-a, Vl - a^, 0) + 6mi + 5^m2 + 0(5^). 
This time, the constraint on the length leads to 



1 = \mf = l+5{—2ami^x + 2\/\ — a^rrii^y) 

+ 6^{-2am2^x + 2Vl - a^ma,^ + .j. + ml y + ml .,) + 0{S^ 
providing the requirements 

(29) Vl^^ 
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Although the 0{S) solution of ruy is no longer zero, these 0{6) terms do not contribute to 
the value of the energy, and the solution to m2,x is required. The constant energy is again 
obtained from the initial conditions of the solution, and is 

We see that the m2,z solution is not needed; it does not contribute to the 0(5^) expansion 
of the energy. 

The expansion of the Hamiltonian system fl20|) leads to the 0{6) set of equations, 

mi,^ = 2(/3^ - (3y)Vl - a? mi,^ 
mi^y = -{-2af3z + K)mi^z 



which are satisfied by 



"T.i,x = cosh cut 

mi , = , smh ut 



where = 4:f3y{(3z — /3y){l — a^). Note these solutions reduce to those in ( l27l) when hx = 
since a = in this case. The 0{S^) expansion of the Hamiltonian system is. 



m2,x = 2(/3^ - - a2 m2,z + '^{Pz - /3y)mi^ymi^z 

m2,y = -{'^Q'l^z + hx)m2,z 

rh2,z = 2Py\/l - a? m2,x + 2(3ymi^xmi^y 

which have the same general solution as the 0{5) equations, but an additional particular 
solution to handle the non- homogeneous part. We only require the solution m2^xi since we 
can then construct the needed m2^y solution from the length constraint f l29|) . Taking another 
derivative of the m2,x equation, and inserting the equation for m2,2, we use the method of 
undetermined coefficients, and obtain the solution 

"^2,x = cosh ut + [{(3^ - /3y) cosh 2ut + sinh 2ut] + - l3y)f. 

Au) 2uj 2 

The initial conditions at 0(5^) are 

m2,x{0) = I + ^{Pz - f3y) 



y/l - m2,y{0) = am2,x{0) 



2(1 - a 



2\ 



Using these, we can construct the expansion of the energy, which after much simplification 
reduces to 



(30) E = Eb- 5^/3y + 0{S'') 

where the saddle point energy, Ei, = /3y + 

The leading order term of the rux component is again used to determine the scaling of the 
period, 

(31) T=- \og^ + 0(1) = - log + 0(1) 
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where = 4/3y(/3;, - /3y){l - hl/4/3^) and Ef, = /3y + This reduces to the period in fl28|) 

when hx = 0. For all values of hx we can use ( l30i) and ( pTi) to describe the energy and period 
near the saddle point in energy, Ei,. 

To finish the approximation of the averaged energy coefficients, we derive the pre-factors 
for their 1/T{E) scalings. Recall that averages are defined as 

1 rTiE) 
if (mm = ^ I fimit))dt 

for trajectories confined to energy E. (Note we will continue to omit the subscripts denoting 
energy branches as we only concern ourselves with branch 1 here.) Having already deter- 
mined the scaling of T as — )■ Ef,, we now turn to the integral f^^^^ f{m(t))dt, which is 
dominated by what happens on the homoclinic orbit with energy Ef,. For the trajectories 
of the components of m{t), we use an approximate trajectory that starts at t = at the 
point on the orbit midway between the two fixed points, where rriz is positive and rriy = 0. 
The trajectories are infinitely long and asymptotically approach the fixed points, therefore 
we approximate 

T(E) poo 

f{m{t))dt^A I f{m{t))dt 



and we take the averages to be approximated by 
(32) (/(m(t))) -J^J^ f{m{t))dt. 

For the simple case when = 0, the exact solution to the Hamiltonian system is 



m^it) = \l — - — ^sechz/ot 



(33) ''TT'yit) = tanhz/ot 



where uq = 2^J f3y{(3z — (3y). Note, this also satisfies E = (3y = fiym^ + /S^m^ and + + 
ml = 1 due to the identity tanh^ x + sech^x = 1. Using ( 132|) to compute averages, we have 



/3z log 1/(5 



^""'^ /3. log 1/5 
{ml) = 1 - (ml) - (ml) 



log 1/6 
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and the average energy coefficients are 



16 
'log 1/5 



^ log 1/5 

For non-zero h^, an exact solution is unknown, but the components rrix and rriz are well 
approximated by sech functions, with appropriate values to match the exact solution at t = 
and as t ^ ±oo. These are 



(35) 



where 



2/3, 

mz ~ ylzsechi/t 



+ AxSechut 



2/3, 



+ 



If3 . - f3y 



\ 



hi 



2/3. 



/3.-/3, 
2Mz 



sechi/t 



hx 



sechz/t 



1/ 



\ 



Ai(3y-hl/APz)i(3z-(3y) + K 



(3z-(3y 



was found by matching the second derivative at t = to the second derivative found from the 
Hamiltonian system. The coefficient was found by noting that wi,(0) = and then solving 
Py + hl/{4:Py) = /3,(1 - ml{0)) - hxTTixiO) for A^. Then, the coefficient A^ = ^/l-ml{0). 
These solutions are also consistent with mx{t) — )■ and rrizit) — )■ as t — oo. Furthermore, 
the solutions in fl35|) reduce to the above exact solutions in fl33|) when h^ = 
For computing the averages, we first note that 

{rrix) = — -j^ + Ax (sechz/t) 



y 



2/3. 

- -^Ax (sechz/t) + A^. (sechVt) 

Py 



and 



(m^) = 1 - {ml) - (ml) 

using the constraint that the magnetization vector has unit length. After computing the 
approximate average defined in ( |32l) . we obtain 



{mx) 



— hrr 



27tA, 



2/3, uTiE) 



hi 



+ 



AAl 



4/32 Pyi^TiE) uTiE) 
AAl 



uT{E) 
. hi 



+ 



2/3, 



vT{E) 
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where T{E) ~ \og{Ei, - E) when E^ = (3y + hl/A(3y and u = 2^ {Py - hi/ A(5y){(5, - ^y). 
While the same procedure could be carried out again for energy branch 2, we point out that 
due to symmetry in the system, we can write the averaged energy coefficients on energy 
branch j = 1,2 (denoted by subscripts) as 



(36) 



where 



ME) 

cm 



uT{E) 
4 

vT{E) 



hi 



Ojhjh,j. + Ti^jhjl3y 
4 



h^ 



uT{E) 



h^ 



6, 




/3. '''2(3y(3, 




hx 



+ K{/3,-(3y 



/3. ^^'^ '"'P. 

(Ji = 1 and (72 = —1, and S = (/3, — Py){Pz — hl/4:f3y). The approximation (136|) are in 
excellent agreement with the averages found via numerical integration, as can be seen in 
Fig. m The limit in Eq. (6) in the text is obtained from the ratios of Aj{E) to Bj{E). 



Appendix D. Matching Conditions 

In this section, we derive the matching conditions for the mean ffist passage time equa- 
tion (7) in the text. Matching conditions are required only at the saddle point where the 
energy branches meet [9] because it is a regular boundary point (see [7] for boundary point 
classification); it is accessible from the interior of each energy branch and the interior of 
each energy branch is accessible from it. On the other hand, no additional boundary con- 
ditions are required at the other ends of the energy branches, specifically the minima, as 
these are entrance boundary points; the interior of the energy branches are accessible from 
these points, but the expected passage time from the interior to these points is infinite. This 
coincides with the diffusion of the magnetization vector on the surface of the unit sphere. 
The original SDE for the magnetization vector contains no extra conditions prescribed at 
the single points corresponding to the energy minima and maxima. The classification of 
the boundary points was determined using the scalings of the averaged energy coefficients 
derived in Sec. [Cl and testing the convergence of certain integrals [7]. 

From the matching conditions, we are able to construct the probabilities of the energy 
switching from one branch to another (the matching conditions required to supplement 
Eq. (2) in the text) as well as the pre-factors for the mean ffist passage times in Eqs. (8) 
and (9) in the text describing the probability the system switches to the other lower energy 
branch rather than return to the original one. The derivation is based on the conservation of 
probability flux of the magnetization vector across the homoclinic orbit on the sphere with 
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energy equal to the saddle point energy, Eh. For ease of notation, any function evaluated at 
energy Ef, should be interpreted as a limit as, E ^ E^ from the interior of the energy branch. 

Consider Pj{E, t) to be the probability density for the energy while on branch j, normalized 
so that 



Pi{E,t)dE+ p2iE,t)dE+ p3{E,t)dE+ p^{E,t)dE = l. 

Ea,l JEa,2 JEt JEt 

These density functions are continuous at the saddle point in energy: 

PliEh,t) = p2{Eb,t) = P3iEb,t) = P4(^b,t). 

The functions Aj{m{t)) are also continuous across the homoclinic orbit, therefore, if this 
orbit is approached from the higher or the lower energy branches, we have that 

fTi(Et) t-T'ziEt) t-T-iiEt) t-TiiEh) 

/ Ai{m{t))dt+ / A2{'m{t))dt= / A3{m{t))dt + / Ai{m{t))dt 
Jo Jo Jo Jo 

or equivalently 

(37) T,{Eh)A,{Eh) + T2{Eh)A2{Eh) = n{Eh)A,{Eh) + T,{Eh)A,{Eh). 

This provides the understanding for why the flux of the total probability, Tj{E)pj{E, t), and 
not simply the averaged probability, pj{E,t), is conserved across the homoclinic orbit. 

The forward Kolmogorov equation for the total probability density, written in terms of 
the flux, Jj[], on each branch j = 1, 2, 3, 4, is 

(38) ^T,{E)p,{E,t) = -^JAT,{E)p,{E,t)] 
where 

JAT,{E)p,{E,t)] =[- aA,{E) + ajB,{E) + 2aeC,{E)'\T,{E)p,{E,t) 

-ae—{A,{E)T,{E)p,{E,t)). 

Analogous to Eq. fl37j) . the conservation of probability flux across the homoclinic orbit pro- 
vides the matching condition for Eq. fl55]) : 

(39) Ji\Ti{Eh)pi{Eh,t)] + J2{T2{Eh)p2{Eh,t)] = JAT^{Eh)pz{Eh,t)] + JAT^{Eh)pi{Eh,t)]. 

The differential equation (7) in the main text for the mean exit time, Tj{E), from energy E, 
comes from the backwards Kolmogorov equation; it uses the adjoint operator to the one in 
fl38|) . Therefore, Eq. (7)'s matching condition is the adjoint condition to the conservation of 
probability flux, Eq. fl39|) . After dividing by ae, the matching condition for Eq. (7) in the 
text is 

A^{Et)T^{Ef,)T[{Ef,)+A2iEt)T2iEhK{Eh) = 

AsiEt)n{EhH{Eh) + A,{Eh)T^{EhK{Eh). 

This condition is equivalent to the condition stated by Fredlein and Wetzell [9]. 

From the exit time matching condition, (HOl) . we derive the probabilities for the energy 
to switch branches in order to complete the stochastic differential equation (2) in the text 
describing the evolutions of the energy, as well as determine the pre-factor for the mean 
switching times between meta-stable states appearing in Eqs. (8) and (9) in the text. 

We deflne the notation P{j — )■ k) to be the probability the energy switches from energy 
branch j to energy branch k at the saddle point, Ej,. In general, this probability is derived 
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from the coefficients of the matching condition (140|) by breaking the integral within the coeffi- 
cients Aj{Ei,) into the parts which lead to each of the other energy branches; these fractional 
parts out of the whole integral yield the probabilities P{j — )■ k). Further simplifications are 
made by taking advantage of the symmetry of this particular magnetic system. 
In general, the probability, P{j k), to switch from branch j to branch k is 

p, . ^ . , Jj^-^"'' Mni(t))h{m(t))dt 

where lk{m{t)) = 1 if m{t) is closer to orbits in branch k than any of the other branches 
besides the one in which it resides, and otherwise. Immediately from Fig. 1(a) in the text, 
we see that 

(41a) P(l ^ 2) = P(2 ^ 1) = P(3 -> 4) = P(4 ^ 3) = 

since lfc(m(t)) = 1 at only two individual points (at the green dots). Exploiting the sym- 
metry about the rrix-my plane, we have that 

(41b) P(l ^ 3) = P(l ^ 4) = 

(41c) P(2 ^ 3) = P(2 ^ 4) = 1, 

(41d) P(4 ^ 1) = P(3 ^ 1), 

and 

(41e) P(4 ^ 2) = P(3 ^ 2). 

Using the above, we can rewrite 

l/2/;^(^^)Ai(m(t))rft _ g. 



(41f) P(3 ^ 1) 



1/2 /J"'^^"^ A^{m{t))dt + 1/2 /o^'^^"^ A2im{t))dt 9i + 92 
where we define the notation 

Qj = / Aj{Tn{t))dt 
Jo 

for j = 1,2. We can take 

^?,^4rf,(/3,2-/3j)- 6,(4/3^2 + 

coming from Sec. IC.2I with out term 4/i/, since gj only appears as fractions. Similarly to 
fl41fj) . we have that 

(41g) P(3 ^ 2) 



9i +92 

All together, the conditions ( 14T|) provide the switching probabilities for the stochastic energy 
equation (2) in the text. 

The mean ffist passage time calculation requires the probability the energy switches from 
branch 1 to 2 or 2 to 1, which we can see from Eq. f l41al) never happens along a direct path. 
Rather, the energy must first switch to one of the two higher energy branches. Conditioning 
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on which intermediate branch the energy switches to, we have that the probabihty the energy 
switches from branch 1 to branch 2 is 

P(switch from 1) = P(l 3)P(3 ^ 2) + P(l 4)P(4 2). 

Using the simphfied probabihties in ( 14T|) . we have that 

r./ u r .^ ^(3^2) P(3^2) g2 
Pfswitch from 1) = ^ H ^ = — — — . 

Similarly, the switching from energy branch 2 back to 1 is 

Pfswitch from 2) = — — — . 

91 + 92 

To match the notation in the text, we define the switching probability from branch j to be 

(42) P(switch from j) = — 

7i + 72 

where 71 = g2 and 72 = f?!- The probabilities in f l42|) are precisely the pre-factors in Eqs. (8) 
and (9) in the text for the mean switching times. 

Appendix E. Mean First Passage Time 

In this section, we derive the mean first passage time equations (8) and (9) in the text. 
Rather than solve Eq. (7) in the text, it is simpler to find the transition time from starting 
point X on energy branch j = 1,2 to the saddle point Eb, then account for the probability 
to transition to the other branch, rather than return to the same well. We therefore find the 
solution, Tj{x), of 

(43) [ - aAj{x) + ajBj{x) + 2aeCj{x)] rj(x) + aeAj{x)Tj{x) = -1 

with absorbing boundary condition rj{Eb) = 0, and multiply it by the switching probability 
in Eq. fH2|) . First we consider the solution valid for any temperature, then consider the limit 
of vanishing temperature. 

The exact solution to f j43l) requires a second boundary condition. As x ^ Ea we know 
from Sec. IC. II that Aj{x) — and Bj{x) — > 0, which leaves the condition 

2aeC{Eay{Ea) = -1. 
Using integrating factors, we integrate fH3|) twice and obtain 



(44) r,(x) = ( ^ ,^ \ , ^ , + -I,{y) ] e(^-^«+^^^(^))A"G^fe)t/2/ 

^ ^ ^ 71 + 72 A \2at{l3y + (3,±K) ae '^'^ ' 

where 



and where 



The expression for 7^ was described in the previous section; it is 

g,^Ad,{Pl-Pl)-h{APl + hl) 
g,^MM-Py)-h{^Py + hl) 

where bj and dj are defined after Eq. ( I36i) . 
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For vanishing temperature (e — )■ 0), rather than approximate fH4|) directly, we notice that 
the solution has a boundary layer where the coefficients Aj{x) and Bj{x) go to zero: both 
near x = Eaj, the minimum (-Ea,i = —hx and i?a,2 = hx), and x = Eb, the saddle point in 
energy. We match the solution coming out of the boundary layer near Ei, to determine the 
leading order expression for the mean ffist passage time from the meta-stable point. This 
meta-stable point is either Eq > Eaj for values of aj when a stable limit cycle exits on branch 
j or the minimum value, Eaj. For simplicity of notation, we will drop the subscript j and 
only consider j = 1. The solution for j = 2 is derived similarly. 

First, we consider the boundary layer near Eb, and rescale the energy by e defining = {x — 
Eb)/e so that ^ — )• — oo leaves the boundary layer. The rescaled equation for g{^) = r{x{^)) 
is 

[ - aAixiO) + ajBixiO) + 2aeC(x(0)] -/(O + = -1 

which to leading order reduces to 

ajB{x{Oy 



-1 + 



a A{x{^)) 

with boundary condition g{0) =0. We then have that 



g'iO = cexp 
and integrating again yields 

9(0 = c / exp 



a 



y 



^B{Eb + ez] 
A{Eb + ez) 

y B{Eb + ez] 



dz 



■dz 



dy 



a Jq A{Eb + ez) 

where we have used the boundary condition g{0) =0. By ffist expanding the integral in the 
exponent in term of e, 

BiEb + ez) 



where 



fy BjEb + ez) ^ roi^^,., a 
L AiE;T7zj'' = l J + Oie)dz^-y, 

A{x) - PI) - h{APl + hi) 
- a J — ^ . 

hhx + TlVhPy (l - ^ j 



\ = a lim 



Bix) 



is the critical spin-polarized current required to induce switching (obtained for the asymptotic 
expansions in Sec. IC.2|) . with h and d given by hi and di after Eq. (!36l) . we have that 



9(0 



Al-aj/\)y^y 



^ A 



,(l-aj/A)C 



1) 



and therefore 

(45) 



1 _ Si 

^ A 



,{l~aj/\){Et,-x)/e 



to leading order. We are left to determine the constant c. As we leave the boundary layer. 



,(l-aj/A)y^^ 



— C 

1 _ oj ' 
^ A 



and we see the solution becomes constant. We turn to consider the full solution in the outer 
region away from the boundary layer to match this constant. 
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Returning to Eq. ( H3|) . and dividing by aeA{x), we have 

-1 ^i_a,BM_,?g(£)V'W + r»- 



a A{x) A{x) J aeA{x) 



e 

We rewrite this as 

(46) re*(^)/V(x)l' = ^^e*(^)/^ 
where 

$(x) = $o(a;) + e<l>i(x) 
for some arbitrary point *. After integrating 0461) from xq to Ef, we have 

The constant from fH5l) . enters through T'{Eh) = g'{0)/e = c/e. Combining with the above 
equation we have 

(47) c = er'(xo)e-(*(^^)-*(^°»/^ - - ^ -i-e-^*^^^)'*^^"/^^?/. 

The integral in ( H7|) is dominated by what happens near Xq, and we have two cases, the 
first when xq is the solution to —aA{x) + ajB{x) = in which case A{xq) ^ 0, and the point 
Xq is away from either boundary layer. The second is when xq is the minimum, and 1/A{x) 
must be canceled by the term generated from the integral of C{x)/A{x) in $(a;). In either 
case, we will need the expansion of 

-iMEh) - $o(x)) = E,-x-^ r ^dy 

a Jx A{y) 

in terms of e defined by ^ = (x — Eiy)/e. We then have 

-iME,) - $o(x)) = -6^ - ^[0 + ee(-l)y + Oie')] ~ -(1 - ^)eC 

a A A 

When Xq is the solution to —aA{xo) + ajB{xo) = 0, the mean passage time, r, is approx- 
imately constant at xq and therefore t'{xq) ~ 0. The expansion of — ($i(i?f,) — $i(x)) only 
contributes higher order terms to the exponent, and 



a A'o ^(y) " J-ooM^o) aA{xo){l-f) 

This, together with r'(xo) = 0, gives the constant in Eq. ( H5l) . Combining with the switching 
probability factor, we have Eq. (9) in the text, 

^ ^ ^ 7i+72«A(i?o)(l-^)^ 

where AE = E^- Eq and £"0 solves -aA(Eo) + ajB{Eo) = 0. 

On the other hand when xq = Ea, the expansion of — ($i(-E'b) — $i(x)) includes a large 
term near ii^^. From the scalings worked out in Sec. IC.lt we know 2C{x)/A{x) ~ l/{x — Ea), 
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which produces a large term, log(x — Ea), in the expansion of — (<l>i(i?b) — $i(x)). Together 
with the scahng A{x) ~ 2(/3j, + + hx){x — Ea) near Ea-, we have 

_g-($(i?,)-<i.fe))A^^ ^ _i / e 



aJE, My) aJE^ 2{l3y + 13, + K){x - Ea) 

1 e-(i-^)« e 1 



For the term in Wl\ involving T'{Ea), we return to Eq. fH3|) . where for x ^ e we have 

2ae(:7(x)r'(x) = -1 

to leading order and therefore 

r'{Ea) 



2aeCiEa) 
We then have 

t'(E )e-(*(^'')-*(^»))/^ ^ lim ^-(^o(En)--S>o(x))/e+\of^(x-Ea) 

" x^Ea 2aeC{x) 

to leading order in the exponent, which goes to zero due to the log(x — Ea) term. Thus, the 
T'{Ea) term does not contribute to the solution. Combining with the switching probability 
factor, we have Eq. (8) in the text. 



7i+72a2(/3, + /3, + /i,)(l-f)2 



(49) r{Ea) ~ :r^-:7T7ir-ir-i-^r^^.^' 

A 

where Ai? = Ei, — Ea- We compare this approximation to the exit time obtained by numer- 
ically solving Eq. fj43|) in Fig. |5l 
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Figure 5. Comparison of the mean first passage time on energy branch 1 
obtained by numerically solving Eq. f l43|) (solid lines) to the approximation, 
Eq. din]) (green dashed lines), plotted as a function of spin-polarized current 
strength, aj (top) and as a function of inverse temperature, e"^ (bottom), 
(top) Between the two thin black lines, the meta-stable starting point is no 
longer the minimum. The constant values = 0.0296 and e = 0.01 have been 
used, (bottom) We have taken, hx = and the spin-polarized currents for the 
lines from top to bottom are aj/a = —0.2, 0.0, 0.2 and 0.4. In both plots, the 
values /3y = 0.06, (3z = 0.5, and a = 0.003 have been used. 
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